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Abstract 

I study the three parameters bipartite quantum Gaussian state called 
squeezed asymmetric thermal state, calculate Gaussian entanglement of 
formation analytically and the up bound of relative entropy of entangle- 
ment, compare them with coherent information of the state. Based on the 
result obtained, one can determine the relative entropy of entanglement 
of the state with infinitive squeezing. 

1 Introduction 

Quantum entanglement is one of the most important phenomenon in quantum 
theory. It exhibit the nature of nonlocal correlation between quantum systems, 
and plays an essential role in various fields of quantum information processing, 
such as quantum computation, quantum communication, quantum cryptogra- 
phy, quantum teleportationpp,and closely related to quantum channel[2]- After 
the first experiments |2j on quantum teleportation using two-mode squeezed 
states @] jS] , a significant amount of work has been devoted to develop a quan- 
tum information theory of continuous variable systems. So far, most of the 
theoretical work has focused on the entanglement properties of the quantum 
states involved in all these experiments, the so-called Gaussian states. The 
first problem arisen is that if a given quantum Gaussian state is entangled, the 
problem of qualifying entanglement has been solved in the general bipartite set- 
ting But the efficiency of entanglement manipulation protocols used in 
practical quantum information processing critically depends on the quality of 
the entanglement that one can generate. It is therefore essential to be able to 
quantify the amount of entanglement in systems with continuous variables es- 
pecially Gaussian states. Several measures have been proposed to quantify the 
amount of entanglement 8 . Among which are the entanglement of formation, 
relative entropy of entanglement and distillation entanglement. Entanglement 
of formation is one of the most important entanglement measures. For symmet- 
ric Guassian state, this entanglement measure can be carried out analytically 
UJ. For general bipartite Gaussian state, another entanglement measure called 
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Gaussian entanglement of formation was introduced^!]]. Clearly it is the up 
bound of the entanglement of formation. Relative entropy of entanglement 11 
[T2*] on the other hand measures the distance between the state under consider- 
ation to the closest separable state. It has the advantage that separable states 
obviously correspond to zero distance. For pure states of bipartite systems it 
reduces to the von Neumann entropy of the reduced state of either subsystem. 
For mixed bipartite states it is usually difficult to be calculated, except for some 
specific states. The distance of the Gaussian state to the set of separable Gaus- 
sian states measured by the relative entropy was considered ^3] ■ Although one 
has yet no proof that there does not exist a non-Gaussian separable state which 
is closer to the Gaussian state under consideration than the closest separable 
Gaussian state, one has good reason to think that it is a good bound of entan- 
glement measure. However, the expression derived ^2] is still not easy for a 
direct numeric calculation. Recently the bound on the relative entropy of en- 
tanglement was calculated for two modes squeezed symmetric thermal state fH], 
In this paper I will consider the more practical used state of two mode squeezed 
non-symmetric themal state, which is to be distributed between two parties by 
means of a lossy optical fiber in asymmetric settings for an initial two mode 
squeezed vacuum state [T5]. I will concentrate on its Gaussian entanglement of 
formation and its up bound of relative entropy of entanglement. According to 
the hypothesis of hashing inequality (TS|, all this entanglement measures should 
be lower bounded by coherent information. I will at last compare the results 
with coherent information of the state. 



2 Squeezed non- symmetric thermal state 

Gaussian state is completely specified by its mean and its covariance matrix, 
where the mean can be dropped by local unitary operation so that is irrelevant 
for entanglement problems. Let us consider a quantum Gaussian state p of two 
system A and B acting on a Hilbert space L 2 (JR.) ® L? (JR.) . (In quantum optical 
term, it is called two modes.) It is convenient to describe Gaussian state density 
operator p by its characteristic function, (e.g. J7j)x ( z ) = Tr (pV (z)) .Here z = 
(.Qa 1 9Bi Pa i Pb ) G 5R 4 is a real vector and V (z) = exp[-z J2k=A,B (<?fcA fc + p k P k )} 
is Weyl operator (displacement operator) with X k and P k are operators of sys- 
tem A and B respectively, satisfying canonical commutation relations . A char- 
acteristic function \ uniquely defines a state p x . Gaussian state is state whose 
X is Gaussian function of z: x( z ) — exp[ir/ T z — jZ T jz), where 7 is a real 
symmetric matrix, called correlation matrix. 

Any Gaussian state of two modes can be transformed into what we called 
the standard form, using local unitary operation only 0. The corresponding 
characteristic function has displacement 77 ~ and the correlation matrix has 
the simple form of four parameters. It was proved that (see ^3 an d reference 
therein, also see [7] and reference therein): For arbitrary real symmetric matrix 
7 there exists linear transformation (symplectic transformation) S p : z — > S p z , 
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preserving canonical commutation relations, then 



7 -> 7o = S p jSp = diag (j A , Jb,Ja, 1b) ■ 



(1) 



This is to say that by proper symplectic transformation, Gaussian state can be 
transformed into a direct product of thermal states p - 
mode thermal state density operators are of the forms 



Po = Po®Po- The one 



Po = 



(1 



OO 

m=0 



(2) 



where v., 



N u +1 



The symplectic transformation 

^u{s p )pu+ (Sp) = Po canBi, 



A r AI = 5(7M-l) ; (^ = AS)- 
5* p induces a unitary operator U (S p ) so that p 
and p = U + (S p ) poll (Sp). For po is a direct product of thermal states. It 
represents the randomness side of state p. While the quantum correlation side 
of p should be caused by the two mode unitary operators U (S p ). 

In this paper I will consider one of the special kind of Gaussian states which 
can be generated from the thermal states po = p$ Pq with a simple form of 



Sp — 



(si 



coshr 
— sinh r 



— sinh r 
coshr 



cosh r sinh r 
sinh r cosh r 



(3) 



Then 7 = S~ ^0 (Sj) \ The Sp 1 induced unitary operator U (Sp 1 ) [H| i 
just the two-mode squeezed operator S2 (r) 



is 



1S2 (r) = exp (a^a^ tanhr) exp [— (a\(iA + a^as + l) In coshr] 

• exp (— aACtB tanhr) . (4) 

This kind of Gaussian states will be called squeezed non-symmetric thermal 
states. And in the rest of this paper, I will use p to specify the states. 

P= S 2 (r) po(v A ,v B ) S} (r), (5) 
Denote A = tanhr, then the inseparability criterion reads 

A > y/VAVB- (6) 

For the convenience of further application, one can express p in coherent 
state representation as 

\ / 2 2 2 2 \ 

(a A ,a B \p\f3A,PB) = C exp[-- [\a A \ + \a B \ + \Pa\ + \Pb\ )] 

■ exp[r (a* A a* B + f3 A (3 B ) + w A a* A /3 A + u B a* B f3 B ] (7) 
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, n (l-v A )(l-v B )(l-X 2 ) \(1-v a vb) v A (l-X 2 ) 

where Cn = , A t — ~r \i , — i — tt, wb = 

Z^l With coherent state representation the reduced state can be eas- 
ily obtained by integration. The reduced density matrix p A — Tr B p and p B = 
TtaP turn out to be one mode thermal states p^ = (l — vjf) J2m=o( v l?) m l m ) mm ( m 
(p; = A, B) with parameters v A d and v r B respectively, 

rd V A (1 - VB) + A 2 (1 - V A ) rd VB (1 ~ Va) + A 2 (1 - V B ) , s 

VA 1-v b + \2(1-v a ) ' Vb 1-v a + \2(1-v b ) ' {0) 

3 Gaussian entanglement of formation 

Entanglement of formation is defined as an infimum 



E F (p) = mi \j2PkE(* k ) 



= ^p fc |^)(^] (9) 



over all (possibly continuous) convex decompositions of the state into pure states 
with respective entanglement being von Neumann entropy of the reduced state. 
By its definition calculating Ep is a highly non-trivial optimization problem, 
which becomes numerically intractable very rapidly if we increase the dimen- 
sions of the Hilbert spaces. Remarkably, there exist analytical expressions for 
two-qubit systems as well as for highly symmetric states. Recently, Ep was 
calculated for symmetric Gaussian states of two modes [0]. But for general two 
mode Gaussian states, it is still not easy if not impossible to carry out the Ep. 
For these reasons the Gaussian entanglement of formation (GEoF) Eq is intro- 
duced to quantify the entanglement of bipartite Gaussian states by taking 
the infimum in only over decompositions into pure Gaussian states. 

For any two- mode Gaussian state with correlation matrix 7 = j q © 7 P , It 
is proved JOj that the Gaussian entanglement of formation is given by the 
entanglement of the least entangled pure state with r ) pur e=X © X~ x which is 
such that 

det (X - j q ) = det (X- 1 - j p ) = 0. (10) 

As a part of the correlation matrix of the entangled pure state, the symmetric 
2x2 matrix X can always be written in the form of 



X 



yx cosh r g y sinh r g 

y sinh r g y x ~ 1 cosh r g 



(11) 



then the pure bipartite state with correlation matrix "f pure can be constructed by 
successively applying two local unitary operations to the two mode squeezed vac- 
uum state, where the two local unitary operations are S p i — diag{y/y, y/y, \J 
and S p 2 = diag{^Jx, V x _1 , Vx -1 , \/x} respectively in simplectic form. 
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Denote n = (N x + N 2 + 1) cosh2r + Ni - N 2 , m = (N x + N 2 + 1) cosh2r - 
Ni + N 2 , k = (N\ + N 2 + 1) sinh 2r. One has 





n 


k 




n —k 


lq = 


k 


m 


1p = 


—k m 



Eq. (|10|l then can be written as 

(n/x + rax) cosh r g — 2k sinh r g — [(nm — k 2 )/y + y] = 0, (12) 
(nx + m/x) coshr g — 2fcsinhr 9 — [(nm — k 2 )y + l/y] = 0. (13) 

So that one has (n + m) (1/x + x) coshr g — 4fcsinhr 9 — (nm — k 2 + l)(l/y + 
y) = 0. Clearly r g is a monotonical increase function of (1/y + y), and when 
y = 1, the minimal value of r g will be achieved. One can subsequently obtains 
x = 1 by substracting Eg. 1121) from Eg. ll3fl . The Gaussian entanglement of 
formation Eg will be 

E G = .g(sinh 2 r 9 ). (14) 
where g(x) = (x + 1) log(x + 1) — a; log a; is the bosonic entropy function, and 

r-=r-ilT,l±^. (15) 

The second term at the right hand side can be written of as Tq, with ^Jvavb — 
tanhro, representing the noise side of the bipartite state. While r represents 
quantum correlation side of the state. The difference of the two gives the GEoF 
squeezing parameter. The compaison of the GEoFs for different ratio of average 
photon numbers Ng/N^ is displayed in Fig. 1. 



4 Up bound for relative entropy of entangle- 
ment 

The relative entropy of entanglement for bipartite quantum state a is defined 
by El: 

E r (a) ee mmS{o \\a) (16) 

where D is the set of all disentangled states, and S (a \\a ) = Tr {a (log a — log 5)} 
is the relative entropy of a with respect to a. Consider the relative entropy of 
entanglement of squeezed thermal state p, if one chooses a subset of D which 
contain all Gaussian separable state, or more specifically all separable squeezed 
thermal state to substitute the district D to carry out the minimum of rela- 
tive entropy, clearly such minimums are local minimums. They can be utilized 
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as up bound E ur (p) of global minimum, the relative entropy of entanglement 
E r (p). Denote the district of all separable squeezed thermal states as Dst, 
then Dst C D, and E ur (p) < E r (p), generally speaking the identity can not 
be achieved. Even for relative entropy of pure Gaussian state, the minimum 
is achieved by non-gaussian separable state Clearly the above way of ob- 
taining up bound by shrinking the district of minimization can be applied to 
other Gaussian state or any other state. In order to obtain the up bound of 
the relative entropy of entanglement of squeezed thermal state p, Let us first 
consider the relative entropy of p with respect to seperable squeezed thermal 
state p(va,vb) — S2 (r) Po (va,v b ) S£ (*")) with A = tanhr < \Jv A v B . The von 
Neumann entropy of state p is |17j 

S(p) = -Trplogp = g(N A )+g(N B ) (17) 

with A,, = j^-, (p = A,B). And one gets log p = S 2 (r) (logp {va,vb))S£ (r) 
by the unitary of S 2 (r). Put pa(v) in explicit operator form, 

Po {v) = {l-v A ){l-v B )rf} aA tf* aB . (18) 
Then the second part of relative entropy will be 

Trplogp = log (1 - v A ) + log(l - v B ) + Tr{S 2 (r) p (v A , v B ) S+ (r) 

^2 (r) {\ogv A a\aA + logiy B a^a B ) S£ (r)} (19) 

By utilizing ^2 (r) ciaS^ (r) = a a coshr — sinhr and with the property that 
operator can be cycled under the trace, after some algebra one obtains 

Trplogp = log(l - v A ) + log(l -v B ) 

+ [- — — — cosh 2 (r — r) H sinh 2 (r — r)] logVA 

l—VA 1 - V B 

+[ - VB — cosh 2 (r — r) H — sinh 2 (r — r)l logws (20) 

l—V B I - VA 

where Tr[(a\a~^ + a,Aa B )po( v A,v B )] = is applied. Let us first find out the 
maximum point of Eq. (|20|l by partial differentiation with respect to r and 
va, v B regardless the fact that A = tanhr < \Jv A v B , it will be at the point 
(r,VA,v B ) — (t,va,v b ). Meanwhile it is noticed that there is no other maxi- 
mum point. Then let us add the condition of A < Vvavb, the maximum should 
be achieved at the edge of Dst, that is A = ^/vav b . After one of the pa- 
rameter is determined, the remain problem is to seek out the maximum with 
respect to va,v b . If the maximum is achieved at (va,v b ) = (y* A ,v B ), and de- 
note y/v A v B = tanh r* , The up bound of relative entropy of entanglement for 
squeezed thermal state will be 
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Eur(p) = -g (N A )-g(N B )- log (l-v* A )+^g(l-v* B ) 

• cosh 2 (r — r*) H sinh 2 (r — r*)] logw^ 



1 - 1 - «b 



cosh (r — r*) H sinh (r — r*)l logUo (21) 

1 - V B l-VA 

5 Hashing inequality and comparison of the 
bounds with coherent information 

One of the most concerned problem in quantum information is the quantum 
capacities, it is shown that the rate of quantum information transformation 
is bounded by the maximal attainable rate of coherent information. If there 
exists hashing inequality, that is, for any bipartite state the one-way distillablc 
entanglement is no less than coherent information, then one obtains Shannon- 
like formulas for the capacities |16| . The coherent information of bipartite state 
a with reductions a a and us is defined as |19j 

I» (a) = S (a„) - S (a) , (22) 

for S (cr^) — S (cr) > and 1^ (a) = otherwise. The hypothetical hashing 
inequality mentioned above is 

D-, (a) > I" {a) (23) 

where (a) is forward classical communication aided distillable entanglement 
of the state. We know that relative entropy of entanglement E r (er) is no less 
than two way distillation of entanglement (cr) [201 > the later is no less than 
one way distillation of entanglement D_> (ct) . Combining with hashing inequal- 
ity one has 

E ur (cr) > E r (cr) > (d) > (a) > I B (cr) . (24) 

One the other hand, it is proved that relative entropy is lower bounded by 
coherent information [21], so E r (cr) > I B (cr) is always true irrespective of the 
hypothesis of hashing inequality. 

Consider the squeezed non-symmetric thermal state p, the reduced states 
p A = Tr B p and p B = TtaP turn out to be not the same. The entropies of its 
reduced states are g(N A d ) and g(N B d ) respectively, where N^ d = v rd /(l — v rd ) 
is the average particle number of the reduced state. This leads to two different 
coherent information I A (p) and I B (p). They are 

I» (p) = max{g(JV; d ) - g(N A ) - g(N B ),0} (25) 
The numerical results are shown in the Fig. 2. 
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6 Conclusions 



A subset of quantum Gaussian mixed state is given and scrutinizingly investi- 
gated. The so called squeezed non-symmetric thermal state has several merits: 
It can be described by three parameters, one for quantum correlation, the other 
two for randomness; The separability of the state readily follows from compari- 
son of the first parameter and the geometric average of the last two parameters, 
if the randomness is stronger, the state is separable, otherwise it is entangled; 
All the three parameters have range from zero to one. It is the simplest Gaus- 
sian mixed state other than the two parameters squeezed thermal state. The 
A, va, vb description of state presented in this paper enable the calculation more 
comparable. With the local unitary operation method, I obtain the GEoF an- 
alytically and write it in a simple style with clearly physical meaning. By the 
nemeric results I find that: As the squeezing parameter tends towards infini- 
tive, the up bound of relative entropy of entanglement tend to coincide with the 
coherent information, so that one can determine the relative entropy of entan- 
glement, and if the hashing inequality is right, the distillation entanglement can 
be determined at infinitive squeezing. 
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Figure 2: Comparison of the GEoF, up bound relative entropy of entanglement 
of with coherent information. With A = 0.99, and Nb/Na = 0.5, dot-dash for 
E ur , solid for I A , and dash line for GEoF. 
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